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ABSTRACT 

We present an analytical model for jets in Fanaroff & Riley Class I (FRI) radio 
galaxies, in which an initially laminar, relativistic flow is surrounded by a shear layer. 
We apply the appropriate conservation laws to constrain the jet parameters, starting 
the model where the radio emission is observed to brighten abruptly. We assume that 
the laminar flow fills the jet there and that pressure balance with the surroundings is 
maintained from that point outwards. Entrainment continuously injects new material 
into the jet and forms a shear layer, which contains material from both the environment 
and the laminar core. The shear layer expands rapidly with distance until finally the 
core disappears, and all of the material is mixed into the shear layer. Beyond this point, 
the shear layer expands in a cone and decelerates smoothly. We apply our model to the 
well-observed FRI source 3C31 and show that there is a self-consistent solution. We 
derive the jet power, together with the variations of mass flux and and entrainment 
rate with distance from the nucleus. The predicted variation of bulk velocity with 
distance in the outer parts of the jets is in good agreement with model fits to VLA 
observations. Our prediction for the shape of the laminar core can be tested with 
higher-resolution imaging. 
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1 INTRODUCTION 

Extragalactic ra dio sources were divided into two morpho- 
logical classes bv lFanaroff fc Rilevl |l973). FRI sources have 
an edge-darkened structure, whereas FRII sources are edge- 
brightened with prominent outer hot-spots. This classifica- 
tion has proved to be extremely robust: the division be- 
tween the classes dep ends primarily on radio luminosity 
l|Fanaroff fc Riievl 1974 ). with FRII sources being more pow- 
erful, but also on the stellar luminosity of the host galaxy 
l|Ledlow fc Owenlll99rj ). There are significant differences be- 
tween the structures of the jets in the two classes: those in 
FRI sources often flare close to the nucleus and have large 
opening angles, whereas their equivalents in FRII source s 
are highly collimated out to the hot-spots (jBridlel 1 1984 ) . 
There is good evidence that FR I jets are initially relativis- 
tic, but decelerate on kiloparsec scales, wh ereas FRII j ets 
remain relativistic until they terminate (e.g. lLaingl [i993h 

The process of deceleration in FRI jets appears to 
be complex, and may involve a transition to turbu- 
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lent flow. In addition, the sources have a wide range 
of morphologies, ranging from well-defined lobes simi- 
lar to those in FRII sources to ex tended plumes or 
tails (|Parma. de Ruiter fc Fantil Il996l h For these rea- 
sons, attempts to construct global models of the evo- 
lution of FRI sources, linking observable quantities 
such as linear size and radio luminosity, have been 
less straightf orward t han t he equivalents for FRII 
sources (e.g. Scheuer 1974 : Kaiser fc Alexander! Il997l : 



iKaiser. Dennett- Thorpe fc Alexander! Il997l h which assume 
that the jet flows are essentially laminar. Part of the mo- 
tivation for the present study is therefore to construct a 
simple model of FRI jets for use as input to global mod- 
els. We consider twin- jet sourc es, which make up at least 



one half of the FR I population ijParma et al 



19961) . exclud- 



ing wide-angle tail and fat -double sources ( Owen fc Laind 
Il989l : [Owen fc Whitelll99lh . whose jet properties differ sig- 
nificantly. 

Over the last few years, detailed modelling of deep 
VLA observations of jets in five FRI sources has al- 
lowed us to quantify their geometries, velocity distribu- 
tions, magnetic fields and emissivity distributions in three 
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di mensions. We re f er in detail to the analysis of 3C 31 
by lLaing fc Bridle! (|2002al . hereafter LB02a) ; observations 
and models of a further four sources have subsequently 
been published (Canvin fc Laingl |2004| ; ICanvin et ail 12003 ; 
lLaing et al.l 12006 ). A consistent picture of FRI jet deceler- 
ation on kiloparsec scales has emerged from these studies. 
The flow velocities are j3 = v/c m 0.8 - 0.9 where the jets 
first brighten abruptly, typically at ~1 kpc from the nucleus. 
The jets flare and then recollimate, decelerating rapidly to 
speeds of /3 « 0.1 - 0.4. The best-fitting transverse velocity 
profiles appear to be approximately self-similar. At least in 
the 4/5 cases where the jets appear to be propagating in con- 
tact with the interstellar medium of the host galaxy rather 
than inside radio lobes, they are roughly 30% faster on-axis 
than at their edges. Nevertheless, an evolution of the veloc- 
ity profiles with distance from the nucleus is not excluded. 
In particular, the transverse velocity variations are poorly 
constrained where the jets first brighten abruptly and a top- 
hat profile would also be consistent with the observations in 
these regions of all five sources. 

In order to decelerate, a jet must entrain matter, either 
from stars within its volume l|Phinnevl 1 19831 : iKomissarovl 
1 1994! ) or by ingestion of the surr ounding mat e rial at its 
bound ary, as originally sugge sted bv lBaanl (|l980h , iDe Yound 
(19811) and iBegelmanl (|l982T l. In the latter case, the trans- 
verse velocity profile almost inevitably evolves with distance 
from the nucleus. 

X-ray observations can be used to infer the temper- 
ature, density and pressure profiles of the hot gas asso- 
ciated with the host galaxies of FRI radio galaxies (e.g. 
Hardcastlc ct al. l l2002l : IWorrall et alh 003: Har dcastle et al.l 



20051 ). Together with the velocity distributions derived from 



modelling of the radio e mission, these can be used in a 
conservation-law analysis fBickncll 1994, hereafter B94) to 
derive jet energy fluxes and the variations of mass flux, pres- 
sure, internal den sity and entrainment rate with distance 
from the nucleus (|Laing fc Bridle! l2002bl . hereafter LB02b). 
Such an analysis is quasi-one-dimensional and therefore 
adopts values for the flow variables (in particular the ve- 
locity) averaged across the jet cross-section. This is reason- 
able if the velocity profiles have restricted ranges and do 
not evolve significantly with distance down the jets, as is 
consistent with the observations of 3C31 (LB02b). If FRI 
jets are in pressure equilibrium with their surroundings af- 
ter they recollimate, this analysis requires that a significant 
overpressure drives the initial flaring. 

An alternative approach, which would also be consistent 
with the observations, is to postulate that the transverse ve- 
locity profiles evolve significantly as the jets interact with 
the external medium. The first approximation is then to as- 
sume pressure equilibrium between the jet and its surround- 
ings and to take explicit account of the interaction between 
the jets and their surroundings using a simple mixing-layer 
model. This is the subject of the present paper. The key as- 
sumption is that there is a turbulent mixing layer between 
the jet and its environment, produced by the interaction 
of the two components. The mixing layer grows both into 
the jet and into the environment, and the initially laminar 
jet eventually becomes fully turbulent. As in th e quasi-one- 
dimensional analysis of lLaing fc Bridie] (|2002bh . we use the 
relativistic formulation of the laws of conservation of mass, 
momentum and energy given by B94. 



We describe the geometry of the jet-layer model in Sec- 
tion [2] The relativistic conservation laws are introduced in 
Section We derive and discuss the solutions for our model 
in Section [4] In Section [5] we apply our model to obser- 
vations of 3C31. We discuss the effects of varying model 
parameters in Section |S] and summarize our conclusions in 
Section [7] 



2 STRUCTURE OF AN FRI JET 

The basic structure of an FR I jet in our model is shown in 
Figure [T] Following the definition given by LB02a, we divide 
the jet into flaring and outer regionsQ Close to the nucleus 
in the flaring region, the outer isophotes have small, but 
increasing opening angles. Further out, they spread rapidly 
and then recollimate. In the outer region, the expansion is 
conical. The radio emission close to the base of the flaring 
region is usually faint and it is always possible to identify a 
distance from the nucleus where the jet brightens abruptly. 
We refer to this location as the brightening poinQ. 

We assume pressure equilibrium with the surround- 
ings at all distances from the nucleus and adopt the sim- 
plest possible p rescri ption for velocity variations following 
ICanto fc Ragal (|l99ll ). Wherever possible, we approximate 
the velocity of a component of the flow by its spatially av- 
eraged value. We postulate that the flow close to the axis 
of the flaring region is laminar, with a constant relativis- 
tic bulk velocity Vj and that this occupies the full width of 
the jet at the brightening point, where interaction with the 
external medium becomes significant for the first time. As 
a result of entrainment of external material, a slower shear 
layer forms between the laminar jet and the environment. 
The bulk velocity of the shear layer should vary continuously 
in the radial (r) direction from Vj at its boundary with the 
laminar core to at its outer edge. The models derived by 
LB02a show, however, that the ratio of edge to centre veloc- 
ity for the synchrotron-emitting material is ~0.7 throughout 
the flaring and outer regions, so the velocity range within 
the mixing layer is fairly narrow. Thus we assume that the 
steady-state flow in this layer has a constant bulk veloc- 
ity v s < Vj. Material from both the environment and the 
laminar jet is continuously injected into the shear layer, the 
latter component supplying energy and momentum as well 
as mass. Integrated across the jet, the fraction of slower ma- 
terial then increases with distance from the nucleus; this 
would be interpreted as deceleration of the entire flow in fits 
to observations with poor transverse resolution. 

The laminar jet in the centre eventually vanishes, so no 
more energy or momentum can be injected into the shear 



1 LB02a postulated the existence of an additional conical inner 
region in the faint inner jets of 3 C31, but observation s of the 
better-resolved source NGC 315 by ICanvin et ail d2005h are in- 
consistent with a constant expansion rate in the corresponding 
part of the brighter jet. A continuously increasing expansion rate 
is required in NGC 315 and is equally consistent with the obser- 
vations of 3C 31 and other sources. A two-zone model is adequate 
to describe the geometry in all cases. 

2 This is also a change of terminology from LB02a, who refer 
to the flaring point, and is done to emphasise that the location 
marks a change in emissivity profile, not in geometry 
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Figure 1. A sketch of the principal features of our jet model (not to scale). For comparison with later figures, the brightening point in 
3C31 is 1.1 kpc from the nucleus and the transition between the flaring and outer regions is at 3.5 kpc (see Fig. [3^,) 



layer from the inside. Motivated by the analysis of 3C 31 
(LB02a), we assume that this transition occurs precisely at 
the end of the flaring region. This may not be general: mod- 
elling of other sources suggests that the bulk of the jet de- 
celeration occurejri_the first part of the flaring region (e.g. 
NGC315; ICanvin et al]|2005h . We assume that the bound- 
ary of the shear layer in the outer region expands smoothly 
and more slowly as the environmental pressure decreases. 
Entrainment from the environment into the shear layer can 
still happen in the outer region, but this requires that the 
velocity be allowed to vary along the jet (Section 13. 2p . We 
assume that there are no transverse velocity gradients. 

The following convention is adopted throughout this 
paper: we use subscript for quantities at the brightening 
point; 1 for quantities at the end of the flaring region; j, I 
and e for all quantities related to the laminar jet, shear layer 
and environment, respectively. Detailed descriptions of the 
parameters are given in Fig. [1] and Table [T] 



3 RELATIVISTIC CONSERVATION LAWS 

We model the structure of FRI jets using relativistic fluid 
mechanics, applying the laws of conservation of mass, mo- 
mentum and energy in the forms given by B94. As in that 
reference, we use the relativistic enthalpy lu = pc 2 + e + p 
and the ratio Si = pc 2 /(e + p) of rest-mass energy to non- 
relativistic enthalpy. Here, e is the internal energy den- 
sity and p is the pressure. In the laminar jet, we expect 



j « P 



-1/4 



(B94). For 3C31, the external pressure drops 
by a factor ~ 4 from the brightening point until the end of 



the flaring region. The approximation that SEj is constant is 
therefore reasonable, and we adopt it in what follows. 

For an ideal gas, e = f/(r — 1), so SS can be written as: 



1 pc 2 



1 mc 



(1) 



where F is the adiabatic index, fh is the average particle 
mass and ks is the Boltzmann constant. is therefore 

a measure of the temperature. We make the approximation 
that the external medium around the jet is isothermal, so 
is constant. There is e vidence for a temperatu re gradi- 
ent on the relevant scales l)Hardcastle et afl 120021 ). but the 
isothermal approximation has a very small effect on our re- 
sults since the energy entrained from the external medium 
is negligible (B94, LB02b) and 3£ e > 1 (Section 15X5) . 



3.1 Conservation laws for the flaring region 

The main difference between our work and that of B94 and 
LB02b is that we divide the flaring region into two parts: the 
laminar jet and the shear layer. Thus our conservation equa- 
tions include distinct terms associated with each of these 
components. 



3.1.1 Conservation of rest mass 

We use the following notations: r j 5 r s are the radii of the 
laminar jet and the shear layer respectively, p is the proper 
density, v is the bulk velocity, /3 = v/c and y = (1 — /3 2 ) _1//2 
is the bulk Lorentz factor. The rest mass of the material 
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passing through the total jet cross section A(x) = nr s (x) 2 
per unit time is equal to the rest mass of the material enter- 
ing through the cross section plus the total entrained mass 
from the environment. We express the mass fluxes M, in the 
laminar jet and the shear layer at distance x separately by: 



Mj(x) = ijPj{x)vjii;ri{x) 2 
M s {x) = r y 3 ps(x)v a n [r s (x) 2 



From equation (9) of B94,we have 
JjPjfiVjirro + / Pe{x')f(x')dx' = Mj{x) + M s (x) 



(2) 
(3) 



(4) 



The first term on the left of equation Q is the rest mass of 
the material entering through cross section per unit time. 
The second term on the left is the entrained mass flux. The 
terms on the right represent the rest masses of the material 
passing through the cross sections of the laminar jet and 
the shear layer per unit time at distance x. We assume that 
the laminar jet continuously supplies energy and momentum 
to the shear layer in such a way that /3j and /3 S remain 
constant throughout the flaring region. The integral term 
gf(x) — f*^ p e (x')f(x')dx / is the mass entrainment function, 
which was given in the form gi(x) = f s pu e nt ■ ndS by B94 
(n is the normal direction of the unit surface dS). f(x) is 
therefore a function that expresses the combination of the 
perpendicular entrainment velocity and the shape of the jet 
boundary. The function gf(x) is a measure of the total mass 
entrained between the nucleus and distance x per unit time. 

We assume that the jet is in pressure equilibrium with 
the external medium throughout the flaring and outer re- 
gions. Thus at fixed x, the pressures in the laminar jet, 
the shear layer and the environment are all equal. Divid- 
ing by p(x) on both sides of equation and defining 
F{(x) = cgt(x)/ [np(x)], we get: 



PO 2 

p(x) 



^£l l3 3 [r s (xf-ri(xf]-F f ( x ). 



(5) 



7 f ^fvrr 2 = Pi(x) + P.(x) + cj>{x), 



(8) 



where cj>(x) = f* g dx' J A (1 - &)ds] is the buoyancy 
term and for pj <C p e in our model, we have 4>(x) = 
fx 7vr e{x') 2 dp. The momentum equation can then be sim- 
plified to: 



PO 2 

p(x) 



1 s — 1 irp(x) 



(9) 



3.1.3 Conservation of energy 

The energy passing through the jet cross section must also 
be conserved. We express the energy flux (or jet power), Q, 
at distance x for the two regions as: 



}j(x) = ^Uj^v-jivrjix) 2 , 



,(x) = j 2 w a (x)v 3 n [r s (x) 2 - rj(x) 2 ] 



(10) 
(11) 



B94 gives the relevant conservation law in his equation (26) 
and we rewrite this as: 



uj e (x')f(x')dx' = Qj{x) + Q s (x). (12) 



As the environment is dominated by the rest mass en- 
ergy, so Sie. is extremely large, and we can approximate 
1 + l/M e w 1 at all position. Thus, /* u e (x')f(x')dx' = 
/* c 2 [1 + 1/M e {x')} f{x')dx' w c 2 g f {x). Dividing both sides 
by cp(x), we get: 



+ i)r, 



r,-i 



Po 2 , ,a 
rtxj r °- r ^ 



' {x) + 1]T Hfr[r s (xf 



r s - 1 



F t (x) 



(13) 



3.1.2 Conservation of momentum 

The momentum flow through the cross section A(x) per unit 
time should be equal to the momentum of the material com- 
ing out of the initial cross section per unit time, modified 
by the effects of buoyancy and differences in pressure be- 
tween the flow and its environment. We express the momen- 
tum flux, P, by: 



Pi(x) = 

Ps(x) = 



nU)Ax] 2 / \ 



Is T^Vs + &PlA X ) 



nrj(x) 2 , 
tt [r s {x) 2 



(6) 
(7) 



where Ap^i(x) =p j (x)—p s (x) and Api, e (x) = p s (x) -p e (x) 
are the pressure differences at distance x. In our pressure- 
matched case, they are all equal to 0. We assume that ma- 
terial is entrained from the environment with a small bulk 
velocity and therefore that it contributes negligible momen- 
tum compared with that of the jet. We rewrite equation (16) 
in B94 for our case as: 



3.2 Conservation laws for the outer region 

For the outer region the conservation equations are similar, 
but without the laminar jet term. Another important dif- 
ference is that the initial cross section is now at the end of 
the flaring region (point 1 in Fig. [TJ. The entrained mass 
and energy now denote the values integrated from point 1 
(distance Xi) up to distance x. Finally, the velocity of the 
layer, /3 S , is a function of distance x. The three equations 
analogous to equations @, IjHJ, and (1121) are then given by 



■yipivmr'i = M s {x) - g a {x), 



li-^vl'Kr'i = Ps{x) + 4>{x) 



•jiujivmr 2 



Qs(x) 



U!e{x')f( X ')dx'. 



(14) 
(15) 

(16) 



The term g (x) — f* p e (x')f(x')dx' is equal to the amount 
of entrained mass per unit time. With the same definitions 
of F[x) and R as given above, these three equations can be 
written in the following ways 
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Table 1. Definitions of key parameters and functions. Columns 4 and 6 indicate whether the values are assumed a priori, inferred 
from fits of relativistic flow models to radio images ('Radio'), derived from X-ray observations of the surrounding hot gas ('X-ray') or 
calculated. 







Flaring 


region 


Outer 


region 


Name 


Physical meaning 


value 


origin 


value 


origin 


r,- 


adiabatic index of the laminar jet 


4/3, constant 


assumed 






I's 


adiabatic index of the shear layer 


4/3, constant 


assumed 


4/3, constant 


assumed 


r e 


adiabatic index of the environment 


5/3, constant 


assumed 


5/3, constant 


assumed 


ft 


bulk velocity of the laminar jet 


constant 


Radio 






ft 


bulk velocity of the shear layer 


constant 


Radio 


function of x 


calculated 


ft 


the bulk velocity at the beginning point of the outer region 






constant 


Radio 


Stj 


ratio of rest mass energy to non-relativistic enthalpy for laminar jet 


constant 


calculated 






Ms 


ratio of rest mass energy to non-relativistic enthalpy for shear layer 


function of x 


calculated 


function of x 


calculated 


Mi 


the value of Ms on the cross section 1 






constant 


calculated 


P 


external pressure on cross section x 


function of x 


X-ray 


function of x 


X-ray 


r i 


the radius of the laminar jet 


function of x 


calculated 






r 3 


the radius of the shear layer 


function of x 


Radio 


function of x 


Radio 


ro 


the jet radius at the brightening point 


constant 


Radio 






n 


the shear layer radius at the beginning point of the outer region 






constant 


Radio 


fit 


entrained mass per time from cross section up to cross section x 


function of x 


calculated 






9o 


entrained mass per time from cross section 1 up to cross section x 






function of x 


calculated 



IVyift Pirj 
Vs - 1 p(x) 
IVri/?iPiri 

r s - l P {x) 



r s7l 2 /3! p x r\ 



m a {x) 



T a - 1 



'js(x)f3s(x)r a (x)' 2 



F (x) 



r s m s (x) + i ,_^ a f ^ 2 _ ,_ N a 



1 p(x) 



r s - 1 Mi + i 

1 0(x) 

Ml + 1 TTp(x) ' 

r s <% s {x) + i 
r s - 1 Mi + i 



-7s (a;) /3s(x) r s (x) + 



7s (x) f3 s {x)r 3 (x) 



, (17) 

(18) 
F {x) 

Mi + r 

(19) 



10" 



4 SOLUTIONS 

In this section, we will solve equations ([5}, (O, (|13|l for the 
flaring region, and equations (|17p . Q18jl. (119[l for the outer 
region in terms of quantities which can be inferred either 
from fits of relativistic flow models to radio images (jet and 
layer velocities in the flaring region, together with the radius 
of the layer in both regions) or from X-ray observations of 
the surrounding hot gas (external density, temperature and 
pressure). We can then derive the shape of the laminar jet, 
rj(x), the variation of velocity with distance in the outer 
region, [3 s (x), the values of Si in the various regions, the 
entrainment function and the velocity of entrainment. 

We assume that the laminar jet has a relativistic equa- 
tion of state with Fj — 4/3; the environment has T e = 5/3. 
The shear layer contains mixed material but the energy den- 
sity must still be dominated by relativistic particles (B94), 
and we therefore take F s =4/3. 



4.1 Solutions for the flaring region 

From equations ((5]l and (|13[) . we have: 



13 
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Figure 2. The external density and pressure profiles for 3C 31. 
The solid lines are derived from the double-beta-model fit to the 
number density and pressure [equations lf33l and {J34j] while the 
dashed lines are power-law approximations with indices of a x = 
1.5 and «2 = 1.1, as described in the text. 
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r a (x) 2 



(&j + 1)7, 



p(x) ' 



'^ s f3 3 {M a {x)-[M 3 (x) + l] 7s } 



(20) 



At the same time, from equation (|9]), we have: 



r a (x) 2 - rj(x) 2 



Mj + l) 7 f$ 



p(x)r a (x) 2 -Por'g 



p(x) ' 



(21) 

Thus, we can express 3£ a as a function of f3j, j3 a and 
the buoyancy term, <j>{ x ), which can be calculated from the 
pressure profile and the shape of the jet r s (x): 

C{x) + B{x) la 



,<%s(x) = 

where 
A = , 



D(x) - A la (3 a 



',-(■«, + 1)7,, 
B(x) = (R 3 + l) 7j /? ; 



- 3 



(22) 

(23) 
(24) 

(25) 
(26) 

Also, from equations ((5JI and (|21|l . we can express the shape 
of the laminar jet and the entrainment function by: 



C(x) = A 



~1sPs + 



3 FjJjPj p(x)irr a (x) 2 - p irr% 

r a - 1 4>{x) 



r s 7 a /3 a p(x)irr a (x) 2 — ponr, 
D{x)=B{x){l- ls ). 



i _ Por 
p(x) 



(27) 



gt(x) - 



N P0 r 






1 p(x)r s (x) 2 


_Q,(x) Qs(x)\ 


, \0 3 t(x) /3 s k(x)] 

[Qs(x) Qj(x) J 





(28) 

where k(x) = [rj(x)/r a (x)] 2 and t(x) = 1 — k(x) are the 
fractions of jet and shear layer, respectively, at distance x. 
Although the expressions for Ps(x) and Pj(x) contain rj(x) 
[equation (|27[) ]. P s (x) /t(x) and Pj(x)/n(x) are functions 
only of observable parameters together with ££ a (x) and Mj. 
By applying the boundary condition rj(xi) — 0, we can de- 
rive and then solve for & a (x) from equation (|22[) given 
the shape of the outer boundary of the shear layer, r a (x). Fi- 
nally, we can determine the shape of the laminar jet bound- 
ary, and the function F(x), which can then be used to cal- 
culate the entrainment function. 



4.2 Solutions for the outer region 

In the outer region there is no laminar jet to supply energy 
to the shear layer but matter continues to be entrained from 
the environment. Thus both (3 a and M s are expected to be 
functions of x. We solve the equations numerically, using the 
following steps. Equations and (|19|l give: 



.(*) = 



ggg[* l(71 -l) +7l ]- 7 .(x) 



7s (x) - 1 
while equation (|18[) gives: 

r s - i A - -Kp{x)(t>(x) 



m s {x) 



r s p s {x)/\m s {x) + i] 



(29) 



(30) 



Again, £% a (x) occurs on the right-hand sides of equation 
and (O, but M a {x)/M s {x) and P s {x)/[M s {x) + 1] are just 
functions of (3 a (x) and other observable parameters. Com- 
bining these two equations, we can solve numerically for the 
value of /3 a (x): the shape of the boundary, r a (x), is con- 
strained from observations, so the only unknown parame- 
ters are /3 a (x), which in turn determines ^f a (x). Then, with 
the known value of /3 a (x), we can express the entrainment 
function as follows. 



9o(x) = 



M 1 7i(*+D-[Sffg + 1 ]^) 



Is (x) - 1 



(31) 



Observations show that the radius of the shear layer in the 
outer region r a (x) increases linearly with x. We use this ob- 
served variation as the input function and predict the dis- 
tributions of /3 a (x), & a {x) and g (x) 



4.3 Summary of the solutions 

To get solutions for both the flaring region and the outer 
region, we adopt the shape function r a (x) from model- fitting 
to radio images, together with the velocities f3 a and j3j for 
the flaring region. We also adopt the pressure profiles from 
X-ray observations. This leaves three functions which need 
to be evaluated at each distance x: rj(x), F(x) and ,'% a (x) 
for the flaring region, and /3 a (x), £% a (x) and Fix) for the 
outer region. 

The three equations from the conservation laws thus 
form a closed system. The input and derived parameters are 
listed in Table [T] 



5 APPLICATION TO 3C 31 

Having established a system of equations which describe 
the structure and kinematics of an FRI jet, we now com- 
pare the results with observational data and models for the 
well-observed source 3C31. Geometrical (projection factor 
and radius) and velocity information are inferred from the 
relativistic-flow models of LB02a. Fits to the density, tem- 
perature and pr essure of the hot gas su rrounding the jets 
are as given by lHardcastle et all (|2002r i and used in the 
quasi-one-dimensional conservation- law analysis of LB02b. 
As in these references, we adopt a concordance cosmology 
with Hubble constant, H = 70kms _1 Mpc -1 , Q A = 0.7 
and Qm = 0.3. At the redshift of the host galaxy of 3C 31, 
z = 0.0169, this gives a scale of 0.344 kpc arcsec -1 . 



5.1 Inferences from observation 

The parameters defining the edge of the shear layer pro- 
jected on the sky are determined by fitting to the total- 
intensity distribution. The angle to the line of sight required 
to correct for projection (52° for 3C31) is derived from the 
relativistic-flow model. In LB02a, the shape of the shear 
layer in the flaring region is described by the polynomial 
r a (x) = a + bx + cx 2 + dx 3 with ro = 0.125 kpc at 1.1 kpc 
and ri = 0.815 kpc at 3.5 kpc. The shear layer initially ex- 
pands slowly, then goes through a phase of faster expansion 
before recollimating at the end of the flaring region. In the 
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Figure 3. Results from our model for the flaring region of 3C 31. (a) Geometry. The outer edge of the flow and the boundary between 
the laminar core and shear layer are shown, (b) Mass flux at distance x. The full and dashed lines indicate the total mass flux and the 
contribution from entrainment, respectively, (c) Profile of M s (x). (d) The entrainment velocity perpendicular to the outer boundary at 
distance x. 



outer region, the shear layer expands conically with an in- 
trinsic half-angle of 13.1°. At the beginning of the flaring 
region, we assume that there is no shear layer, so we use 
the on-axis bulk velocity inferred by LB02a to characterize 
the jet, Vj — 0.77c. We suppose that the shear layer makes 
up essentially all of the flow at the end of the flaring region. 
LB02a infer a variation of velocity across the flow from 0.37c 
- 0.55c at this distance so we take a representative value of 

Vs = 0.45c. 

lHardcastle et all (|2002l ) have estimated the external 
density and pressure profiles for 3C 31 from X-ray obser- 
vations. The density profile is given by: 



p e {x) = m p n e (x)/xH, 



(32) 



where m p is the mass of a proton, \h = 0.74 is the abun- 
dance of hydrogen by mass and n e (x) is the proton number 
density of the environment given by: 

re. (a) = n c (l + x 2 /x 2 )- 3 ^ /2 + n g (l+ x 2 /x 2 g )~ 3 ^ /2 . (33) 

The numerical values of the parameters are: n c — 1.8 x 
10 5 m" 3 , n g = 1.9 x 10 3 m- 3 , /3 C = 0.73, (3 g = 0.38, 
x c = 1.2 kpc, x g = 52 kpc. The temperatures estimated by 
lHardcastle et alj (|2002l ) range from 4.9 x 10 6 K to 1.7 x 10 7 K, 



corresponding to ffl e — 5x 10 5 to 1.5 x 10 5 . Thus the approxi- 
mation l + 1/.^e ~ 1 (Sec tion 3.1.3) is valid to hi g h acc uracy. 
The pressure is given bv lBirkinshaw fc Worraul (|l993l ): 



p(x) = kT{x)n e {x)/(nxH), 



(34) 



where p = 0.6 is the mass per particle. For simplicity, we 
approximate the pressure and density distributions using 
power- law forms: 



p e {x) = p e , { — )' 
Xo 

p(x) =p { — )~ a2 

Xo 



(35) 
(36) 



where xq is the position of the brightening point. p e 



2.16 x 10" 



;m and po = 1-93 x 10 Pa are the den- 



sity and pressure at xo, respectively. The values Qi = 1.5 
and Q2 = 1.1 give good approximations to the profiles, and 
we adopt them in the following calculation. The correspond- 
ing density and the pressure profi les are compared with 
those from lHardcastle et al.l (|2002| ) in Figure [2] Although 
we use an isothermal approximation in the development of 
our model (Section [3}, the assumed pressure profile includes 
the effects of the temperature gradient. 
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Figure 4. Results from our model for the outer region of 3C31. (a) Profile of bulk velocity s (x). (b) Mass flux M(x) at distance x. 
The full and dashed lines indicate the total mass flux and the contribution from entrainment, respectively, (c) Profile of M s (x). (d) The 
entrainment velocity perpendicular to the outer boundary as a function of distance, x. The jagged shape of the profile is a numerical 
artefact, but the overall shape is correct. 



5.2 Results from the model 

5.2.1 Flaring region 

With the parameters given in Section 15.11 we obtain Blj = 
13.4 in the flaring region. The profiles of M s (x) and the 
total mass flux passing through a given cross section, M, 
are plotted in Figure In the same figure, we also plot iw, 
the normal component of the entrainment velocity across 
the surface of the jet. This is related to the entrainment 
function by v cnt = (1/ p e )dg/ds. 

The model predicts that the laminar jet initially ex- 
pands at the beginning of the flaring region and then starts 
to collapse ~1.7kpc away from the brightening point. Mean- 
while, the value of M a [x) drops a little at the beginning of 
the flaring region and then reaches an asymptotic value of 
ss6.7. The initial decrease of M s {x) occurs because the small 
amount of entrained material at the beginning of the flar- 
ing region can easily be heated by the the laminar jet. The 
functional forms of M a (x) and v Bn t(x) are constrained by 
the parameters inferred for 3C 31 and may differ in other 
sources. For example, if the shear layer initially expands 
faster, R 3 (x) will be higher and v ent lower throughout the 
flaring region. 



5.2.2 Outer region 

In the outer region, our model predicts that the bulk velocity 
P a should decrease smoothly with x. @ s = 0.45 at 3.5 kpc, 
where it is normalized to the mean value of the distribu- 
tion derived by LB02a, decreasing to 0.22 at 12 kpc. This 
is reasonably consistent with the velocity range derived by 
LB02a (P = 0.15 - 0.22 at the same distance). The value of 
Ms increases with x in our solution, reflecting the increasing 
dominance of the mass by entrained material. We plot M s {x) 
and p s (x) together with profiles of mass flux and velocity in 
Figure [4] 



5.2.3 Estimate of jet power 

Using the calculated and observed parameters given above, 
we can estimate the power of the jets in 3C 31. The relevant 
para meter for comparis on with estimates by other methods 
fe.g. lBirzan et al.ll2008r i is $ (LB02b), the energy flux of the 
jet with the rest-mass contribution subtracted. $ = Q — Mc 2 
in the notation of the present paper. Applying equation (|10[1 
at the brightening point, we get values of Q — 3.4 x 10 37 W 
and $ = 1.6 x 10 37 W. 3C 31 is a fairly powerful FRI source, 
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Figure 5. (a) The entrainment function, g(x), from our model 
(full line) compared with the estimate from stellar mass loss 
within the jet, g a (x) (dotted), (b) As in panel (a), but for the 
entrainment per unit length of the jet, dg/dx. 



with a monochromatic luminosity of 10 24 ' 5 W at 1.4 GHz, 
approximately a factor of 10 below the FRI/FRII dividing 
line plotted by iLedlow fc QwerJ d 19961) given the absolute 
magnitude of the host galaxy (|Owen fc Laing||l989r ) . A total 
power of $ = 1.6 x 10 3r W for the twin jets of 3C 31 is well 
within the range derived from observations of cavities in the 
X-ray gas surrounding oth er radio galaxies o f comparable 
monochromatic luminosity l)Birzan et al.ll2008h . 



5.2.4 Mass input from stellar mass loss 

It has been argued that the deceleration in the flaring region 
could be caused by the entrainme nt of stellar wind m aterial 
from stars located inside the jet ( Komissarovl Il99~4l ) . In or- 
der to test this idea, we adopt the estimate of mass input 
from LB02b, who u sed a deprojection of R-band surface pho- 
tometry for 3C 31 |o wen fc LaingHl989l ). together with the 
same assumptions o n conversion betwe en st ellar luminosity 
and m ass loss as in iKomissarovl fl994) and iBowman et al.l 
(1996). The corresponding entrainment per unit length (the 
derivative of the entrainment function defined above) can be 
written as 



1.0 F 
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X, 

Figure 6. The relation between Mj and /3j for the flaring region. 
The values of ro, ri, p(x) and f3 s are fixed at the values determined 
for 3C 31. The plus sign indicates the value of Sij for 3C 31. 
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Figure 7. The relation between and j3 s for the flaring region. 
The values of ro, r\, p(x) and /3j are set to the values determined 
for 3C 31. The plus sign indicates the value of Mj for 3C 31 



dg s /dx = 2.4 x 10 28 7rr s (:r) 2 a;~ 2 - 65 kgkpc" 1 yr~\ (37) 

where r a (x) and x are in units of kpc. In Fig. [5] we compare 
the entrainment function from our model and its deriva- 
tive with those estimated for stellar mass loss. At the be- 
ginning of the flaring region, the stellar mass input rate is 
remarkably close to that required, given the crudity of the 
assumptions. At larger distances, however, it falls well below 
the level required to decelerate the jet. In the outer region, 
the entrainment rate per unit length required by our model 
continues to increase, whereas that from stellar mass loss 
decreases. Thus, although stellar mass loss may be impor- 
tant in initiating the jet deceleration at the start of the flar- 
ing region, boundary-layer entrainment, as described by our 
model, is clearly required on larger scales. Mass input dis- 
tributed throughout the jet volume, as would be expected 
from stellar mass loss, is a potential complication to our 
analysis. 
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Figure 8. The relation between and ct2 for the flaring region. 
The values of ro , r\ , /3j and /3 S are fixed at the values determined 
for 3C 31. The plus sign indicates the value of Mj for 3C 31. 



Table 2. Comparison between derived parameters for 3C31 in 
this paper and LB02b. Following B94 and LB02b, we quote the 
rclativistic Mach number, ^# = fvv/-f Cs Cs, where c g is the sound 
speed and 7 Cs = [1 — (cs/c) 2 ] -1 / 2 . 



Quantity 


This paper 


LB02b 


Fnerev flux (TO 37 W) 


1.6 


1.1 


(excluding rest mass) 






Initial momentum flux 


7.7 


3.7 


(10 28 kgrns" 2 ) 






Density at brightening point 


12 


2.5 


(l(T 27 kgm- 3 ) 






Mass flux at brightening point 


6.2 


1.0 


(lO^kgyr- 1 ) 






Mass flux at 12kpc 


47 


32 


(lO^kgyr- 1 ) 






Pressure at brightening point 


1.9 


15 


(1CT 11 Pa) 






M at brightening point 


13.4 (jet) 


0,1 




7.7 (layer) 


0,1 


Mach number at brightening point 


7.7 (jet) 


1.5 




2.5 (layer) 


1,5 



5.3 Comparison with LB02b 

It is of interest to compare the results of the present model 
with the conservation-law analysis of LB02b. The treatments 
are very similar in many respects, both relying on quasi-one- 
dimensional approximations and using conservation of mass, 
momentum and energy in a realistic external environment. 
The formulation of the conservation laws is identical in the 
two treatments. The principal differences in the assumptions 
are as follows. 

(i) The analysis of LB02b explicitly assumed that there 
are no variations in physical parameters across the jets, as 
in our treatment of the outer region. We split the flaring 
region into laminar jet and shear layer components. 

(ii) The jets in LB02b's analysis are assumed to come into 
approximate pressure equilibrium with their surroundings 
only after they recollimate. This then requires that they are 
over-pressured at the start of the flaring region. In contrast, 
we assume that the jets are everywhere in pressure equilib- 
rium with the external medium. In this picture, the initial 
expansion is caused by transfer of momentum from the lam- 
inar core to the shear layer rather than a pressure- driven 
expansion. 

(iii) The models are constrained in slightly different ways. 
Both specify the radius of the jet as a function of distance 
from the nucleus. In LB02b, the velocity is given everywhere 
and the best average match to pressure equilibrium is found 
for the outer region. In the present model, velocities are 
specified only in the flaring region, but pressure equilibrium 
is enforced along the entire length of the jet. 

(iv) In the solutions preferred by LB02b, momentum flux 
= $/c initially. This is required for the jets to decelerate 
from highly-relativistic velocities on parsec scales, as in uni- 
fied models of BL Lac objects and FR I radio galaxies. It is 
not an explicit constraint in the present models, where the 
momentum flux is relatively higher (corresponding to the 
solutions in section 3.3.6 of LB02b). 

(v) We use power-law, isothermal approximations for the 
external density and pressure distributions, whereas LB02b 



use a double-beta-model with varying temperature. The dif- 
ferences are minor (Fig. [2} ■ 

LB02b discussed the effects of varying the assumptions of 
their analysis. This led to a spread of values around those for 
their reference model which we quote here. Table[2]compares 
values of key parameters for our model jet and that from 
LB02b's reference model at the brightening point and at 
12kpc from the nucleus. 

The energy fluxes of the two model jets are quite similar, 
despite the differences in starting assumptions. In terms of 
the available energy flux <E> (with the rest-mass component 
subtracted, as in Section 15.2.31 and LB02b) , we find $ = 
1.6 x 10 37 W, compared with $ = 1.1 x 10 37 W for LB02b. 
This is because the geometries of the two jets are identical; 
in the outer region their velocities are very similar and they 
are both close to pressure equilibrium with the surroundings. 
The main difference is in the mass flux, which is a factor of 
1.5 times larger at 12kpc from the nucleus in the present 
model. 

There is a larger difference between the initial condi- 
tions for the two models at the brightening point. The model 
jet of LB02b has an initial density roughly 5 times lower than 
that described here, but is also overpressured: its energy den- 
sity is dominated by the internal energy of relativistic par- 
ticle rather than by bulk kinetic energy, as can be seen from 
the differences in the value of & at the brightening point 
(Table [2]). The very low initial density in LB02b's reference 
model is derived from the requirement for FR I jets to be able 
to decelerate from bulk Lorentz factors ~5 on parsec scales. 
If this requirement is relaxed, as in the high-momentum so- 
lutions described in section 3.2.6 of that paper, results closer 
to those in the present paper are obtained. The entrainment 
rate at the beginning of the flaring region in both models 
is very low and could be provided by mass input from stars 
(Section l5.2.4|l . Both models require an additional source of 
mass at larger distances from the nucleus, however. 



Relativistic mixing-layer model for FRI jets 11 



6 EXPLORING PARAMETER SPACE FOR 
THE MODEL 

Our model uses several parameters derived from observa- 
tions of 3C 31 to calculate the key physical properties of 
this object. For other FRI sources, these parameters will be 
inappropriate and in this section, we discuss the effects of 
altering them. 

6.1 Flaring region 

The parameters affecting the solution in the flaring region 
are the value of 3£ e , the polynomial coefficients for the outer 
boundary, the jet and layer velocities and the gradient of 
the external pressure. We have argued that M e , which is 
always very large, cannot affect our solutions significantly. 
The shape of the outer boundary plays an important role in 
determining the buoyancy term and varies from source to 
source. As the shape function has four free parameters, we 
will not discuss this point in detail her^E but we note that 
faster expansion of the shear layer will lead to larger values 
of M s (x) and smaller entrainment velocities. We can vary 
the remaining three parameters, /3j, f3 3 and 02, individually 
to determine their effect on our solutions and we plot them 
against Mj below. The distributions of M s , mass flux and 
v en t are closely related to that of @t$ . 

Given that the laminar jet is assumed to be in pressure 
equilibrium with its surroundings at the brightening point, 
its internal energy is determined. If (3 S and the form of the 
pressure profile are also fixed, then the energy flux minus 
the rest mass term, $ (defined by its value at xi) is also 
unchanged. Since $ is a conserved quantity, this is also true 
for the laminar jet at Xq. A faster jet with the same internal 
energy must therefore have smaller density and Rj (Fig. [HJ . 

Moreover, if we have a faster shear layer at Xi, which 
means that $ is higher, but fij remains constant, then the 
density of the laminar core at xo must increase, since the 
internal energy is fixed there by the pressure balance con- 
dition. Rj therefore increases with f3 s (Fig. 0. The shapes 
of the distributions of gt(x), & s (x) and v ent (x) remain the 
same but their normalizations change if the jet or layer ve- 
locities are varied. For a faster laminar jet or a slower shear 
layer, 3£ s and Vent (x) both become smaller, indicating that 
the shear layer is less dense. 

The value of Mj also depends on the pressure profile, 
quantified here by the exponent oti of a power-law distri- 
bution. If the pressure decreases more slowly with distance, 
then the assumption of pressure equilibrium requires the in- 
ternal energy of the layer to be higher at the end of the 
flaring region, increasing the energy flux. If the velocity of 
the laminar core is fixed at the brightening point, as is its 
internal energy, then we need a denser laminar jet and there- 
fore a higher value of Mj (Fig. [8}. 

6.2 Outer region 

For the outer region, the situation is much simpler. As , /3i 
and )"i are determined by continuity at the boundary with 

3 More recent model s use a two-parameter form for the shape 
of the flaring region llCanvin fc Laing|l20oj ; ICanvin et all 120051; 
lLaing et ai1l2009) . 



the flaring region, the only additional parameters inferred 
from the observations are the half opening angle 9 and the 
power- law exponent of the external pressure profile, 0:2- Two 
factors influence the opening angle: the decrease of external 
pressure and the expansion associated with entrainment. Of 
the two, the latter is more important for 3C31: if we set 
Vent = to remove the entrainment terms, the predicted 
jet opening angle is around 3° (compared with the observed 
value of 13°), suggesting that entrainment dominates the 
expansion. 

Figure [9] shows how the jet properties change as func- 
tions of the exponent of the external density and pressure 
distributions, 02. For a jet with a fixed opening angle, a 
larger value of 02 (a faster decrease of pressure) reduces the 
amount of material entrained from the environment into the 
jet and leads to a slower entrainment velocity. As the buoy- 
ancy force can accelerate the material in the jet, a larger 
value of Q2 can also lead to a slower deceleration in the 
outer region. The outer region cools due to continuous en- 
trainment of thermal matter from the environment into the 
shear layer, so M a increases with distance at a rate depen- 
dent on the entrainment velocity. 

If we keep 0:2 = 1.1 and alter the opening angle, 6, the 
jet properties vary as shown in Figure [lOl We find that when 
the opening angle is small, the jet hardly entrains any mate- 
rial from the environment, and so decelerates more slowly. In 
extreme cases, the jet could even be accelerated slightly by 
the pressure gradient. It is interesting to note that the other 
four sources which have been mod elled in detail al l have 
outer region opening angles < 5° (|Canvin fc Laind |2004| ; 
ICanvin et ail 120051 ; lLaing et af] I2006T ) and show little evi- 
dence for deceleration on these scales. Compared with 3C 31, 
their external environments are significantly less dense and 
it may be that entrainment is relatively less important at 
large distances from the nucleus. 



7 CONCLUSIONS AND FURTHER WORK 

We have constructed an analytical mixing-layer model for 
jets in FRI radio sources that satisfies the relativistic mass, 
momentum and energy conservation laws. FR I jets are ob- 
served to expand rapidly and then recollimate into conical 
outflows, and we divide them into flaring and outer regions 
based on this morphological distinction. We assume that the 
jet is in pressure equilibrium with its surroundings through- 
out both regions and divide the flaring region into two parts: 
a laminar jet with very high bulk velocity, and a slower shear 
layer. We prescribe the shape of the shear layer and the (con- 
stant) velocities of the laminar jet Vj and shear layer v s in 
the flaring region. We can then derive the jet power Q and 
the ratio of rest mass energy to non-relativistic enthalpy for 
the laminar jet, R s . We calculate profiles along the jet of 
the mass flux M{x), the entrainment velocity v e nt(x) and 
the ratio of rest mass energy to non-relativistic enthalpy for 
the shear layer, R 3 (x). Finally, we predict the variation of 
the bulk velocity of the shear layer, v„(x), with distance from 
the nucleus in the outer region and the radius of the laminar 
core r s (x) in the flaring region. 

We have applied the model to the well-observed FRI 
radio source 3C 31, and find self-consistent solutions for the 
jet properties. In the flaring region, we take the shape of 
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Figure 9. The jet properties in the outer region for different values of ce2, the exponent in the external pressure distribution. The solid 
line is the value estimated for 3C31, «2 = 1-1. The dotted line, dash dot line and dashed line are for «2 = 0.5, 02 = 1-5 and a.2 = 2, 
respectively, (a) Velocity profile, /3 s (x). (b) The entrainment function g {x). This is the entrained mass flux between the start of the outer 
region (x = xi) and distance x. (c) Profile of M s (x). (d) The entrainment velocity perpendicular to the shear layer surface. Irregularities 
in the profile are numerical artefacts. 



the shear layer r„(x) and the bulk velocities of Vj — 0.77c 
and v a — 0.45c from fits to VLA observations (LB02a). 
In the outer region, our model predicts that the bulk ve- 
locity should decrease smoothly to 0.22c at 12kpc, which 
is consistent with the values derived by LB02a. The cor- 
responding energy flux is Q = 3.4 x 10 37 W, equivalent to 
$ = 1.6 xlO 37 W if the rest -mass contribution is subtracted. 

We find that Mj = 13.4 and that M a [x) in the shear 
layer decreases from ~7.5 at the beginning of the flaring 
region to 6.7 and then stays almost constant until the jet 
recollimates. In the outer region, 8£ 3 {pi) increases from 6.7 
to 15.7 at 12kpc, indicating that the temperature of the 
material in the outer region is decreasing with distance. The 
velocity of entrainment into the jet varies with distance, but 
has a characteristic value of a few hundred ms~ . 

Our model gives a somewhat larger energy flux for 3C 31 
than that of LB02b, who find <E> = 1.1 x 10 37 W assuming 
that there are no transverse velocity variations in the jets. 
The two models are quite similar in in the outer region, but 
differ more significantly at the start of the flaring region: 
our analysis assumes pressure equilibrium whereas LB02b 
require a significant over-pressure and consequently find a 



lower initial density. Both models require entrainment rates 
which are consistent with estimates of mass input from stars 
at the base of the flaring region, but not at larger distances. 

We plan to apply a slightly generalized version of our 
analysis to the other FRI jets for whic h velocity mod- 
els and adequate X-ray data are available dCanvin fc Laind 
|2004| ; ICanvin et alJliooBI ; iLaing et al.ll2006T ). Complex, non- 
axisymmetric structures are observed at the start of the 
flaring regions of these jets, as they are in 3C 31 (LB02a). 
It is plausible that these are shocks in the supersonic flow 
required in the core, although the detailed morphology 
of the best-resolve d example, NGC 315, suggests otherwise 
|Laing et al.ll2006fl . Our model requires that there should 
be a clear demarcation in velocity between the core and the 
shear layer in FRI jets and predicts the shape of the for- 
mer. This can in principle be tested using the techniques 
developed by LB02a, but existing observations are limited 
by insufficient resolution or sensitivity in regions of rapid 
deceleration close to the nucleusQ EVLA and e-MERLIN 



4 Transverse velocity gradients are clearly detected, but they are 
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Figure 10. The jet properties in the outer region for different values of the opening angle, 9. The solid line is the default value for 3C 31 
with 9 = 13.1°. The dotted line, dash dot line and dashed line arc for 9 = 3.5°, 9 = 8° and 9 = 20° respectively. 



should be able to image the flaring regions in detail and to 
resolve a core/shear-layer structure if one is present. 
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